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Abstract 

There are currently no models readily available that provide nucleon-nucleon spin dependent 
scattering amplitudes at high energies (s > 6 GeV^). This work aims to provide a model for calcu- 
lating these high energy scattering amplitudes. The foundation of the model is Regge theory since 
it allows for a relativistic description and full spin dependence. We present our parameterization of 
the amplitudes, and show comparisons of our solutions to the data set we have collected. Overall 
the model works as intended, and provides an adequate description of the scattering amplitudes. 
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I. INTRODUCTION 



We present a model for calculating elastic, spin-dependent scattering amplitudes for the 
nucleon-nucleon system. While much work has been applied to this topic over the years there 
is no analysis available for both proton-proton and proton-neutron in the mid to high energy 
range, Mandelstam s > 6 GeV^. Thus far the most complete, highest energy, and readily 
available work is the SAID [U [2] analysis which provides the proton-neutron amplitudes to 
s ~ 6 GeV^ and the proton-proton amplitudes up to s ~ 9.8 GeV^. Our goal is to calculate 
the amplitudes at higher energies. In order to accurately describe the nucleon-nucleon system 
at these energies we require a fully relativistic, spin dependent model. Furthermore, due to 
the scarcity of data, particularly in the proton-neutron case, we require a model which will 
provide some confidence in extrapolating the results to higher energies. 

Our primary motivation in building this model is to utilize the amplitudes to describe the 
final state interactions of deuteron electrodisintegration. It has been shown that a complete 
description of the final state interactions is necessary in order to accurately describe this 
process [3H5]. Currently the final state interactions are given by the SAID amplitudes, 
however, the kinematics at Jlab, where experiments have been performed, allow for final 
state nucleons with energies greater than can be described by SAID. 

Our approach is to parameterize the helicity amplitudes in terms of Regge poles or ex- 
changes P-[TU]. Regge theory has had great phenomenological success over the years. The 
theory is fully relativistic and allows for a complete spin-dependent description. Regge 
theory is also ideal for us to use since it readily scales to higher energies. Furthermore, 
Regge theory has been utilized in the past to model proton-neutron scattering at mid-range 
energies with good results [llj. Overall Regge theory provides us a systematic method of 
parameterizing the scattering amplitudes, while meeting all the criteria of the model. 

The fundamental idea of Regge theory is to study the analytic behavior of the ampli- 
tudes, when one allows the angular momentum J to be continuous and complex. While the 
analysis is rigorous [B] for non-relativistic scattering, the relativistic case is based on a series 
of assumptions. However, it is in the relativistic case, by exploiting crossing symmetry, that 
Regge theory is useful as a parameterization method. Performing the Regge analysis in the 
t-channel center of momentum frame, the large s, small t approximation to the amplitudes 
in the s-channel center of momentum frame is obtained [3^9] • 
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While Regge theory excels at high energies, at lower energies it becomes more difficult 
to implement, as more and more Regge exchanges can contribute. Because of this feature, 
however, it naturally lends itself as a method for extrapolating to higher energies, since as 
one increases in energy the low energy poles are suppressed. In our approach we are able 
to fit to the low energy nucleon-nucleon data and extrapolate our results to higher energy 
regions where data is unavailable. 

Section [TT] discusses our method of parameterizing the helicity amplitudes in terms of 
Regge poles. Then in Section ITlT] we present comparisons between the data and our results 



for both a polarized and unpolarized fit solution. We conclude in Section IV with a summary 
and outlook. 



II. THEORETICAL FRAMEWORK 

All observables in the nucleon-nucleon system can be described by five independent he- 
licity amplitudes [T2]. These amplitudes are given as. 
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In order to calculate the amplitudes we look to Regge theory to provide us with a param- 
eterization method. Applying Regge theory to the nucleon-nucleon system presents some 
challenges, primarily due to the inclusion of spin. The Regge analysis should be performed in 
the crossed (t) channel, and the result crossed back to the direct (s) channel. Because of the 
many helicity configurations the crossing relations are complicated. Also, Regge exchanges 
have definite quantum numbers, such as parity P, G-parity G, and isospin /, which must 
be taken into account, and because of the symmetries of the nucleon-nucleon system, any 
non-strange mesonic Regge exchange can contribute. Finally, since nucleons are fermions 
we need to properly take into account Fermi statistics, shown in Fig. [TJ 

Fortunately we can can either avoid, or at least simplify, these complications by relating 
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FIG. 1. Pictorial represention of the helicity amplitudes in terms of the Fermi invariants. P 
represents the various gamma matrices which contribute to this process. 

the Regge exchanges to the Fermi invariants [131 IE! ' 

f = Fi(.,t)l«.l(^)-F/(s,t)(z75)(^)-(n5)(^^ 

+ F^(s,t)7'^«7f + Fi{s,t)i^,YY'\^,^,Y'^ (2) 

where s and t are the Mandelstam variables, / is an isospin label, 1 and 2 correspond to 
the vertices shown in Fig. [1} This immediately benefits us since we get all spin dependence 
"out in the open" , and we can immediately do the Dirac algebra to get the s-channel helicity 
amplitudes, 

TPP^PP ^ J2 {C*. [Ffis,t) + Fl{s,t)] - Cr, [f!;{s,u) + Flis,u)] } (3) 
j 

Tr^^- = J2{C'^^[F^is,t)-F}{s,t)] -2C-F}{s,u)}, (4) 
j 

where i corresponds to the helicity configurations (++; ++), (++; H — ), (H — ; H — ), (++; ), 

(H — ; — h), and j to the different types of Fermi invariants S, V, T, P, A. The matrices C*^- 
and Cj^-, containing all the spin dependence, are obtained from performing the Dirac alge- 
bra, and are given in the appendix. For convenience Mandelstam u is used in the terms 
corresponding to the interchange of the final state particles necessary to account for Fermi 
statistics. Fig. [l|b). 

Our goal, instead of Reggeizeing the helicity amplitudes directly, is to instead parame- 
terize the Fermi invariants in terms of Regge exchanges. This is extremely beneficial since 
the crossing relations become trivial as the Fermi invariants are invariant, and since we have 
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taken care of spin explicitly, we will see that the Regge analysis will reduce to the spinless 
case. 

A. Relation Between Fermi Invariants and Regge Exchanges 

Regge exchanges are found in the crossed channel and have definite quantum numbers 
P, G, and I. I is easily factored out and is taken care of in § and Q, so we simply need 
to retain the label here. Also, we need only concern ourselves with exchanges related to Fig. 
[TJ^a), since Reggeization of Fig. can easily be obtained by interchanging t u in our 
final result. In order to find the Regge contributions to the nucleon-nucleon (A^A^ — )■ A^A^) 
system, we need to analyze the t-channel, nucleon-anti-nucleon (A^A^ — )■ A^A^) amplitudes of 
definite P and G. 

We will focus only on the initial state since there is a simple relation between the initial 
and final states. Our basic state that we will work with is, 

iAn)al3 = Vai-pt,\2)Ufi{pt,\l)\I,Mi), (5) 

where pt is the t-channel center of mass momentum, Ai and A2 are the helicities of particles 
1, and 2, I is the total isospin of the state, Mj is the 3'''^ component of isospin, and we have 
labeled the Dirac indices explicitly with a and /3. The goal is to symmetrize this state in 
terms of parity and G-parity. We begin with parity. In Dirac space the parity operator is 
7o. Parity acting on this state yields, 

P{i'in)aP = -Va{pt, -\2)up{-pt, Mi) (6) 

We can construct a vertex with definite parity then as, 

ii^^U^al = ^ H-Pt, A2)r«n(p-;, Ai) - Pvip-t, -A2)r«w^(-pi, -ao] |/, Mj) (t) 

Defining 7or7o = PrF, where Pr is ±1, we can simplify ([T]) to, 

(Oa/^ri^j = ^ii + PrP)v{-Pt^2)r^Mpt^i)\i,Mj) (8) 

We now move on to G-parity. G-parity is defined as, G = Ce'"'^^^, where 12 is the y rotation 
matrix in isospin space and C is the charge conjugation operator, given in Dirac space as, 
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TABLE I. Symmetries of 7 matrices, 
r / 7^ 7*^ 7 7'^7^ 77^ id s 
Pr + - + - - + - + 
Cr + - - - + + - - 

C = C'foK, where K is the complex conjugation operator, and 

\-icx2 J \' / 

acting on a two particle state, (|5|, yields, 

G'(^/'i„)a,3 = Ua{-Pt, >^2)Vl3{Pt, Ai ) (- 1) V2^Ai^C ^z) (10) 

where r]±\ = (— l)^/^^'^, and r]c is an arbitrary phase which is convenient to define as 
Tjc = (— l)^r7A2^Ai- We can then construct a vertex of definite parity and G-parity as. 



(V'f/)«^r2 = ^(1 + PrPM-pt, A2)r(i)n(pl, Ai) + Gui-p-t, \2)T^%ip-u |/, M,) (11) 
Defining CpF = C^yoKTC^yoK, where Cr is ±1, this simplifies to, 

ii^fr^'UTl^l = ^il + PrP)il + Vx,Vx,CrG)v{-p-'t,\2)T^'^uipt,\i)\I,Mj). (12) 
Pr and Cr for the available couplings are given in Table |Ij Note that we use the decomposition 



The result, (12) can also be used to calculate outgoing states of definite parity and G- 



parity by utilizing the relations, 

{u{p, \)Tv{-p, \2)T = vH-p, \2)r^u^p, Ai) = v{-p, \2)ru{p, Ai). (13) 

With these results we can construct the symmetrized amplitudes, 

^.fAt;A„A. = Sl^'F^is, t) + P^f'FUs, t) + V^^F^{s^ t) + Aj^Fi(., t) + Tl^Fi{s^ t), (14) 

where Aj represents Ai, A2, A'^^, A2, which are the helicities of the initial and final particles, 
and S^l^, Pxf^, V^^'^, A^^^, are obtained from the dirac algebra, and are dependent upon 
specific values of P, G, and helicity configurations. 
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Now, in order to Reggeize, we set up a partial wave expansion, of definite parity and 
G-parity, in tlie t-channel, center of momentum frame. 



(15) 

where Et and 6t are the t-channel center of momentum energy and scattering angle, and the 
fy^x'^ Xi \2^^^) correspond to partial wave coefficients from expanding on to the Wigner-d 
functions. 



We now equate (14) and (15) and select specific P, G and helicity values to isolate and 



solve for each of the Fermi invariants in terms of partial waves. Essentially, there are many 
redundant equations since the Fermi invariants do not depend on spin, giving us the freedom 
to only choose helicity combinations of (++;++) and (H — ; H — ). Once this is done all the 
Wigner-d functions reduce to Legendre polynomials, and the Reggeization procedure reduces 
to the spinless case. Each Fermi invariant trivially crosses back to the s-channel, ([s]) and 
(|4]), and we do not have to deal with any of the complications that are associated with a 
typical Regge analysis of particles with spin. From here we can follow the typical methods 
to Reggeize spinless amplitudes [7H9], the result is. 
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(16) 



where m = .93895 (GeV) is the nucleon mass, and the right-most vector is defined by a sum 
of Regge exchanges. 



RTf{s,t) = C(^,t) ^a±(t)/3r^(t) -1 



2s 
4m2 - t 



(17) 
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where (3(t) and a{t) correspond to the residue and the trajectory of the Regge pole and are 
discussed in the following section, (i^^t) is a cutoff factor also discussed in the following 
section, j is simply the position of the Regge exchange in the vector, and 



Ut) = { 



g-i{7ra{t)/2+<5) _^ 

(18) 



Here we have also introduced an additional phase for each Regge exchange, S, which 
accounts for the fact that we have absorbed all extra t dependence into the residue, including 
the approximation we utilized for ^±(t). Ultimately it provides an extra degree of freedom 
which is convenient when fitting certain Reggeons. Also note that while Reggeons with 



PG = , enter into two different positions in (16), the residues of any contributing poles 

in these positions are not necessarily the same. 

In order to take into account the u-channel exchanges of Fig. [l](b), we can simply sub- 
stitute t — )■ u in (16). We also utilize an additional factor of for type 4 exchanges. 



guaranteeing that amplitude c? = at t = 0, which is required by conservation of angular 
momentum. In addition, we multiply type 5 exchanges with a factor of which we as- 
sume we can factor from ^4(5,^). This is necessary in order to cancel with an additional 
factor of s in the matrix C^j and prevents amplitude e from blowing up at large s. This 
seems to be a general problem with expressing the amplitudes in terms of the Fermi invari- 
ants at large s, and Fa{s, t) should either always be redefined or parameterized to explicitly 
cancel this factor of s in order to avoid this problem. Now that we have the Fermi invariants 
parametrized in terms of Regge exchanges we can plug this result into ^ and ^ for a 
Regge approximation of the s— channel helicity amplitudes. 



B. Residue and Trajectory 

We utilize linear Regge trajectories, a{t) = ao + ait. These are obtained from the meson 
masses available from the Particle Data Group[15]. In addition to the mesonic trajectories, 
we also utilize "effective" trajectories as discussed below. 
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We use three different parametrizations for the residues, 

Pn{t) = (1 - e^*) f3oe^^' (19) 

where Pq, Pi and 7 are fit parameters. We utihze the different types of residues for different 
Regge exchanges, as well as different fit solutions, based on trial and error. 



Equation (17) differs from the usual expression in that we have kept the full expression 
for cos(6't). Generally the Regge limit assumes that cos{6t) » 1 which impies that s » 
Am? — t. In extrapolating from the region where the SAID partial wave analysis has been 
performed to higher s, we are violating this condition in two respects. First, data where s 
is of the same order of magnitude as Am? are included. Second, in the same region there 
are significant data for 4m^ — s < t < 0. So at backward angles t is of the same order of 
magnitude as s. For this reason we keep the exact expresssion for cos(6'(). 

We have chosen to fit the low s data using the same form as the Regge parameterization, 
but in doing so it is necessary to introduce additional "effective" trajectories that are not 
derived from the meson spectrum. 

A practical problem associated with fitting these forms at low s is that the u channel 
contributions necessarily overlap those from the t channel. Fitting to data near 6' = 0°, 
where t = and u = 4m^ — s, and near 6 = 180°, where u = and t = 4m^ — s, can be 
affected substantially by the tail of the crossed channel. This can cause the fitting procedure 
to become very sensitive, if not unstable. As a result we have found it useful to introduce a 
cutoff factor, 

C(s,t)= fl-eKsi^-Oy (20) 



to decouple the t and u channel contributions at the endpoints in order to simplify the fitting 
procedure. This has no effect at large s where the two channels have no significant overlap, 
but is extremely useful for smaller values of s. 



C. Electromagnetic Effects 



In order to properly describe the proton-proton interaction we need to account for electro- 
magnetic effects. We use the full proton vertex (El), with a one photon exchange. The one 
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photon exchange amphtudes are given in the appendix ( E4 ) . In order to account for higher 
order effects we utihze a hehcity-dependent constant and phase. Since the electromagnetic 
contribution is dominated by "no flip" and "single flip" contributions we redefine the one 



follows, 






(21) 


h'EM{s,t) = (5be'HEM{s,t) 


(22) 




(23) 



where /3a, /3b, /?c, ^a-, ^b, and Sc are fit to available polarization and differential cross section 
data. Utilizing this approach allows us to keep the electromagnetic effects under control, 
and smoothly fit between the Coulomb and hadronic regions. 

III. RESULTS 

We present here two solutions to our fits, an unpolarized version as well as a polarized 
version. While one of our primary goals is to include all spin dependence, we also foresee 
applications for which an unpolarized solution will prove useful. 

In comparison with the SAID analysis, there is overlap between the two models at lower 
energies, and we expect the SAID solution to be more precise, i.e. lower x^- The emphasis 
on our fit is the ability to extrapolate to higher energies. Using a Regge model over the 
entire angular region may give less precise results than the SAID parameterization, but it 
does allow this extrapolation. A further test of our results will be when we implement them 
into the electrodisintegration of the deuteron process, as it will allow us to see how much 
the results vary in comparison with the use of the SAID amplitudes through the kinematical 
region of overlap. 

The data set was assembled from the SAID analysis [2], the Durham database [16], the 
Cudell dataset [17j, and the Particle Data Group We also reference here the original 
papers [TBl - 1184] . The dataset that we have collected can currently be obtained by contacting 
the authors. 

All observables were fit simultaneously. The values are given in Table O, and the 



parameter values are given in Tables III and IV for the polarized and unpolarized solutions 



respectively. In order to avoid the largest data set dominating the fit, we implemented 
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weights which were varied in order to keep all observables on the same footing. This was 
especially useful since the proton-neutron data is so limited in comparison to the proton- 
proton data. Because our data set includes differential cross section data from many sources, 
there is a potential problem with normalization. In order to correct for this we fit to the 
shape of the differential cross section data and allow the overall magnitude of the data to 
float by plus or minus 15%. 

TABLE II. values for both the unpolarized and polarized solutions. The unpolarized solution is 
based on 6111 data points and has 132 parameters. The polarized solution is based on 13869 data 
points and has 138 parameters. 
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FIG. 2. Total cross sections for proton-proton and proton-neutron as a function of Mandelstam s. 

The total cross sections for both proton-proton and proton-neutron are presented in Fig. 
[2j As these are calculated at t = 0, the Regge approximation works extremely well. 

In order to constrain the model at large s, we also fit to high energy proton-proton data. 
In Fig. [3] we show differential cross sections through both the Coulomb and dip regions, as 
well as the polarization parameter. These results illustrate the ability of the Regge model 
to scale to higher energies. 

Low energy differential cross sections for proton-proton and proton-neutron are shown 
in Figs. |4] and |5| The model works very well, especially considering that we describe the 
data over the entire angular region, and for relatively low s, well outside of where one would 
typically expect the Regge approximation to be valid. We have also included the unpolarized 
results in these graphs. While the unpolarized results yield a very low there is limited 
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FIG. 3. High energy results for proton-proton differential cross sections in both the Coulomb and 
dip region, as well as high energy polarization. 



data to constrain the fit in the pn case. 

Single polarization or analyzing power are presented for proton-proton, Fig. [6} and for 
proton-neutron Fig. [7j Again the model describes the data well, although more proton- 
neutron data would be useful to constrain the model further. 

Finally we present the double-polarization observables, proton-proton in Fig. |8] and 
proton-neutron in Fig. [9j These were fit with minimal priority, due to the lack of data. For 
each of these observables we roughly describe the data, and in certain cases the model works 
very well. 
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d (deg) 

FIG. 4. Differential cross section for proton-proton as a function of center of mass angle 9. Each 
data set is offset by a factor of two. 

IV. SUMMARY AND OUTLOOK 

We have parameterized the nucleon-nucleon scattering amphtudes in terms of Regge 
exchanges. Relating to the Fermi invariants allows us to calculate all spin dependence 
directly, while ensuring that a Regge exchange with definite quantum numbers contributes 
appropriately to the amplitude. Because of the ability of Regge theory to scale, we are able 
to extrapolate our results to s ~ 20 GeV^ with reasonable confidence. 

Our next step is to utilize these amplitudes to describe the final state interactions of the 
D(e,e'p)n process. We also plan to provide these amplitudes, as well as the data set we have 
collected to the nuclear physics community. Currently the amplitudes and data set can be 
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FIG. 5. Differential cross section for proton- neutron as a function of center of mass angle 9. Each 
data set is offset by a factor of two. 



obtained by contacting the authors. 
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FIG. 6. Polarization for proton-proton as a function of center of mass angle 9. 
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FIG. 7. Polarization for proton-neutron as a function of center of mass angle 6. 
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FIG. 8. Double polarization observables for proton-proton as a function of center of mass angle 0. 
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FIG. 9. Double polarization observables for proton- neutron as a function of center of mass angle 
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Appendix A: Parameters 



Parameter values for polarized and unpolarized results are given in Table |III[ and Table 



IV respectively. 



Appendix B: Amplitudes and Observables 

All observables can be written in terms of the five independent helicity amplitudes 
j;iven in ([T]). We present here the observables relevant to this paper, 

-2m2 



a 



^J s{s — Am?) 

^ = , "\ + 4|6P + |cp + MP + |ep) (B2) 
at 27rs(s — 4m^) ' 

a=i(|a|2 + 4|6|2 + |cp + |rf|2 + |e|2) (B3) 

da _ 
dt 7rs(s-4m2) 

aP = aAN = -'^[b*{a + c + d-e)] (B5) 

aAxx = ^ia*d + c*e) (B6) 

aAzx = -^[b*{a + d-c + e)] (B7) 

aAzz = -l{\a\' + \d\'-\c\'-\e\') (B8) 

aAyy = ^{a*d - c*e) + 2\b\^ (B9) 

aD = ^{a*c-d*e) + 2\b\'^ (BIO) 

a^T = 3^(a*e - d*c) + 2|6|2 (Bll) 



Appendix C: Helicity Spinors 

In the center of momentum frame the helicity spinors are, 

u{±p,\)=N\ ^ ]x±x{p), (CI) 

v{±p,X)=N\ U^a(p), (C2) 
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TABLE V. Two component spinors of (CI) 
xi(p) X_i(p) 



initial state 



final state 



c 


1 1 











/ COS ^ 1 


M 


f - sin ^ \ 






^ cos| / 



where = 
in Table |Vl 
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p is a unit vector in the direction of p, and x±a(p) are given 



Appendix D: Amplitudes to Fermi Invariants 



The helicity dependent matrices which relate the Fermi invariants to the helicity ampli- 
tudes are, 
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Appendix E: Electromagnetic Effects 



We utilize the full proton vertex given as, 



GE{t)-GM{t)tlAm^ GM{t)-GE{t) 
GE = GM/2.79={l-t/.71)-\ 



(El) 
(E2) 
(E3) 



where q^, is the four momentum of the photon. 
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The one photon-exchange contribution to the hehcity amphtudes is then given as, 

I \ 47r/137 
"^^^^''') = 2t(4m2-.)(-4m3 + mt)2^ ^^^^ 

(-SGijCA/m^stn - G\t(?,2m^ + sH + 2m2t(s + t) - 8m^(s + 2t)) (E5) 
- 8G|m^(16m^ + 25^ + 3st + - 4m2(3s + 2t))) (E6) 



47r/137 / stu 

^^^^ = ~ 2mt(-4m2 + t)^ \/ (. - ^rn^f ^ ^^^^ 

((s - M)(4m2G'| + Gl^t) + 2G'£;G'M(16m^ - - 4m^(s + t))) (E8) 

/ X (47r/137)n 

GEM s,t = - ^ w . 3^ .sa X E9 

2t(4m^ — s)(— 4m'^ + mt)"^ 

(8G|m^(u - s) + 8GijG'Mm2st + G'y(-8m'^ + 2mh - st)) (ElO) 

, / X 47r/137 , 

^""^^'^^= (.-4m^)(t-4rnT ' ^''''^ 

(4GbGa/sm + ^15(5 - m) + GM(16m^ + 23^ + 3st + - Am'^{?,s + 2t))) (E12) 

/ N 47r/137 , 

{AGeGmsu + ^13(5 - m) + G^(16m^ + 2s^ + 3st + - 4m2(3s + 2t))). (E14) 
where u = Am? — s — t. 
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